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The Lyman-Qf emission of high-z damped Lyman-ar systems 
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ABSTRACT 

Using a spectral stacking technique we searched for the average Lya emission from high-z 
Damped Lya (DLA) galaxies detected in the Sloan Digital Sky Survey QSO spectra. We used 
a sample of 341 DLAs of mean redshift < z > = 2.86 and log N(H i)>20.62 to place a 3cr 
upper limit of 3.0xl0~ 18 erg s cirT 2 on the Lya flux emitted within ~1.5 arcsec (or 12 kpc) 
from the QSO line of sight. This corresponds to an average Lya luminosity of < 2 x 10 41 
erg s _1 or 0.03 L*(Lya). This limit is deeper than the limit of most surveys for faint Lya 
emitters. The lack of Lya emission in DLAs is consistent with the in situ star formation, 
for a given N(H i), being less efficient than what is seen in local galaxies. Thus, the overall 
DLA population seems to originate from the low luminosity end of the high redshift Lya 
emitting galaxies and/or to be located far away from the star forming regions. The latter may 
well be true since we detect strong O vi absorption in the stacked spectrum, indicating that 
DLAs are associated with a highly ionized phase possibly the relics of galactic winds and/or 
originating from cold accretion flows. We find the contribution of DLA galaxies to the global 
star formation rate density to be comparatively lower than that of Lyman Break Galaxies. 
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1 INTRODUCTION 

Damped Lya (DLA) systems, detected in absorption in the spectra 
of background quasars are characterized by large neutral hydrogen 
column densities (N(H i) > 2 x 10 20 cm -2 ) and repre sent the main 
reservoir of H i at high-z (see Noterdae me et al.l2009l Fig. 14). The 
presence of associated heavy elements, the evolution with redshift 
of the mass density in DLAs, a signature of gas consumption via 
star formation and the detectability of DLAs over a wide range of 
redshift make them the appropriate laboratories for studying the 
cosmological evolution of star formation activity in a luminosity 
unbiased way. 

Though observational studies of DLAs have been pursued 
over 25 years, one of the most important questions that remain 
unanswered yet is the connection between DLAs and star-forming 
galaxies. Measuring the star formation rate (SFR) in DLA-galaxies 
is very important as DLAs could provide substa ntial contributions 
to the global SFR d ensity at high redshifts dWolfe et all [2003b : 
ISrianand et alj2005l) . 

At low and intermediate redshifts (z < 1), galaxy coun- 
terparts of DLAs have been identified i n a number of cases 
dBurbidge et aill 19961 : iLe Brun et alii 19971 ; iRao et al.|[2003h. They 
are u sually low surface brightness dwarf galaxies (see Rao et al.1 
l2003h . On the other hand, searches for the direct emission of the 
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high redshift (i.e z > 1) DLA galaxies have resulted in a num- 
ber of non-detections (jLowenthal et alJI 19951 iBunker et al.|[l999l: 
IColbert & Malkanl 120021 : iKulkarni et al j 12004 IChristensen et all 



2007) and only a few cases have been spectroscopically con- 
firmed through dete c tion of Ly« emissio n (Warren & Moller 1996; 
i Moller et al.1 1 19981: IMoller et al.1 I2002L |2004 iHeinmuller et all 
l2006l : lFvnbo et aljboiCh . The difficulty is mainly attributed to the 
faint nature of the DLA galaxies and their apparent closeness to the 
bright background QSOs. 

One indirect way to address this question is to study the rel- 
ative populations of different rotational levels of molecular hy- 
drogen (H 2 ) together with that of fine- structure levels of Ci and 
C n, The ambient UV flux and therefore the in-situ SFR can thus 
been derived (e. g. Hirashit a & Ferraral 2005 ; ISrianand et aHl2005l : 
iNoterdaeme et al.ll2007h in a few of the ~ 10% of DLAs that show 
detectable amounts of H2 JPetitiean et alj200dl : lLedoux et al.l2003l : 
Noterdaeme et al. 2008). The rotational excitation seen in these H 2 
bearing DLAs is consistent with an amb ient UV radiation fiel d sim- 
ilar to or higher than the Galactic one JSrianand et alj[2005T) . This 
may not be representative of the overall population however. 

Another indirect approach is to derive the SFR through the 
cooling rat e inferred from the C 11* A 1335 absorption lines as sug- 
gested by IWolfe et alj J2003al) . They have estimated an average 
SFR per unit area of 10~ 22 M Q yr 1 kpc~ 2 or 10" 13 M Q yr~' kpc~ 2 
if DLAs arise in cold neutral medium (CNM) or warm neutral 
medium (WNM) respectively. However, appreciable contribution 
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of C n* absorption from the warm ionised medium (WIM) as ob- 
served in the diffuse gas of the Milky Way dLehner et al]2004l) will 
make the inferred SFR based on Cn* absorption lines unreliable. 

Thanks to the availability of a very large number of moder- 
ate resolution QSO spectra in the SDSS data base, it is possible 
to directly detect emission from the galaxy counte rparts within 1.5 
arc sec to the line of sight dNoterdaeme et alj|201dh and/or to mea- 
sure the average SFR in different types of QSO absorpti on systems 
by detecting nebular emission in the stacked sp ectra ( Wi ld et al.l 
l2007l : lNoterdaeme et alJ|201ol ; lMenard et al.ll2009l) . Here, we tried 
to measure the average Lycr emission from DLAs using spectral 
stacking techniques. We describe our DLA sample in Section[2] dis- 
cuss the stacking technique we use and the results in Section[3] We 
discuss the O vi detections in Section[5]and contribution of DLAs to 
global SFR density in Section [4] Our conclusions are presented in 
Section[6] Throughout this paper we assume a standard flat ACDM 
universe with H = 71 kms~' Mpc~\ O m = 0.26, £1 A = 0.74. 



2 DLA SAMPLE 

The sample of DLAs used here mainly comes from an automatic 
se arch for DLAs in the Slo an Digital Sky Survey II, Data Release 
7 dNot erdaeme et al. 2009). This sample contains 1426 absorbers 
at 2.15 < z < 5.2 with logMHi) > 20, out of which 937 sys- 
te ms have log jV(H i) > 20.3. We also use additional DLAs found 
by IProchaska & Wolf el J2009h . In total, there are 914 DLAs with 
redshift measurements based on metal absorption lines. We notice 
that the poor background subtraction leaves some spikes in the lo- 
cations of strong sky lines. To avoid any systematics due to these 
features we have excluded systems located in regions of strong sky 
lines Q. 

In order to establish the presence or absence of any Lyo- emis- 
sion, we need the bottom of the DLA absorption profile to be 
defined by enough pixels with zero flux. We will use only those 
DLAs for which the core (with zero flux) extends over at least three 
FWHM elements. Since the spectral resolution of SDSS spectra is 
about R = 2000, the rest FWHM of an unresolved Lyo- emission is 
~ 0.61 A. The above condition translates to a lower limit on the H i 
column density of log N(H i) = 20.62. 

By visual inspections of the individual spectra, we rejected 
systems showing a double absorption feature in either their Ly/J or 
metal absorption profiles. Indeed some of these systems may not 
be real DLAs as the high Lyo- equivalent width may be the result of 
the blend of several components. We also rejected proximate DLAs 
located within 5000 km s -1 from the QSO redshift and ensured that 
there is no contamination of the sample by broad absorption lines. 
After applying all these conditions we are left with 341 systems 
with log N(H i)>20.62 that forms our sample to be used in the 
stacking exercise. 



3 SPECTRAL STACKING 

The observed SDSS spectra were first shifted to the rest frame 
of the DLA conserving the flux per unit wavelength interval. The 
residual flux level in the core of the DLA absorption profile gen- 
erally shows some non-zero off-set probably related to sky sub- 
traction errors. We correct this off-set using the median flux in the 

1 The wavelength range affected by sky lines are 4040-4060A, 4340-4380 
A, 5170-5230 A, 5450-5475 A and > 5560A. 



Table 1. hya measurement in the stacked spectrum. 

Stacking method F(Lyaf F 3lr (Lya) h L 3cr (Lyo-) 6 

simple mean 0.21±0.13 039 283 

weighted mean 0.26±0.13 0.39 28.3 

weighted mean with clipping 0.13±0.10 0.30 21.8 

" in units of 10~ 17 ergs s cm -2 
* in units of 10 40 ergs s _1 



absorption trough before stacking. While this exercise makes our 
stacked spectrum insensitive to the presence of continuum light 
from the DLA galaxy, it will not affect the detectablity of the Lyo- 
line. 

The flux in each pixel of the combined spectrum is calculated 
by averaging the fluxes of all spectra at this position. Arithmetic 
mean and weighted mean are both used. In the case of weighted 
mean the i-band signal-to-noise ratio (SNR) is used as the weight- 
ing factor for each spectrum. We calculate the uncertainty in each 
pixel of the combined spectrum as the standard deviation around 
the mean value. We obtain a third average spectrum from the 
weighted mean of the values in each pixel after rejecting the 5% 
extreme values on either positive and negative sides. 

The stacked spectra obtained using simple mean, weighted 
mean, and weighted mean with clipping together with associated 
errors are shown in Fig.[JJ We also show the stacked synthetic ab- 
sorption DLA profiles in the corresponding samples as red contin- 
uous curves. The dashed blue curves are the synthetic profiles ob- 
tained using the lower and upper limits on N(H i) obtained in indi- 
vidual systems. To calculate the Lya emission flux, in each stacked 
spectrum, we restrict ourselves to the wavelength range covered by 
the core of synthetic spectrum obtained from the lower limits on 
A'(Hi) (i.e inner dotted profiles in Fig. [TJ. We obtain the limit on 
the hya emission line flux by integrating the observed flux over a 
gaussian function with FWHM of 200 km s -1 . This value corre- 
sponds to the mean velocity widths of low ionization lines in DLAs 
dLedoux et aT]|2006h convolved with the SDSS instrumental broad- 
ening. The observed FWHM of the low ionization lines in our com- 
posite spectrum is consistent with this value. 

The Lyo- flux measurements are summarised in Table[JJand the 
gaussian fits are shown in the insets of Fig.Q] We do see marginal 
Lyo- flux at the expected position at a < 2a level in the case of 
simple and weighted mean. In the case of clipped mean the bottom 
of the hya absorption profile is consistent with no Lyo- emission. 
We estimate the 3<x upper limit (F 3cr (Lya) given in Table [TJ by 
using the gaussian fitting errors. As expected, the 3cr limit from 
the clipped mean is smaller than that from the other two meth- 
ods. As this is the case with less contamination by outliers, we 
use the upper limit from the clipped mean composite spectrum, 
Fio- (hya) < 3. Ox 10~ 18 erg s _I cirT 2 , for all further discussions. For 
the SDSS fibre diameter, 3 arcsec, this limit translates to a limiting 
lcr surface brightness limit of 1.4x10"" ergs s cirT 2 arcsec -2 . 
This is only 1.7 times higher than tha t achieved in the v ery deep 
long slit spectroscopic observations by Rau ch et al] d2008l) . 

The flux limit we have achieved in the stacked spectrum is 
almost one order of magnitude smaller than the flux measured in 
the cas e of direct detections of hya emissi on from DLA galax- 
ies (see lM0ller et alJl2004l : lFvnbo et alj|20ld) a nd the typical flux 
limits reported for individual measurements (see lChristensen et al.l 
2007). In addition the flux limit we have reached is much deeper 
than the one reached in typical blind narrow band searches for Lyo-- 
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Figure 1. Stacked spectrum of DLAs obtained from simple mean, weighted mean, and weighted mean with 5% clipping are shown in, respectively, the left, 
middle and right panels. The red continuous curves are the stacked spectrum of the DLA synthetic absorption profiles. The dashed blue curves are the synthetic 
profiles obtained using the lower and upper limits on N(H I) as measured in each DLA. A gaussian with FWHM = 200 km s~' is fitted in the bottom of the 
absorption trough and is shown by a red continuous curve in the inset. Top panels: profile of C n/11334 line to give an idea of the velocity spread of neutral gas. 



emitters (see lOuchi et al.l2008h and co mparable to the deepe st limit 
achieved in the VVDS-LAE survey bv lCassata et alj d201dh . 

For the mean redshift of our sample, < z >= 2.86, our limiting 
flux corresponds to a Lya- luminosity of 2x 10 41 erg s -1 (see TableQJ 
or ~ 0.03L* (Lya) if we use the Schechter function par ameter L* = 
5.8 x 10 42 erg s -1 from the narrow band Lya search bv lOuchi et al.l 
(2008). Thus it appears that on an average, DLAs at z ~2.8 trace 
the faint end luminosity of Lyo- emitters. 



4 STAR FORMATION 

We estimate the average SFR in DLAs assuming that Lyo- photons 
mainly originate from H n regions around massive stars and case 
B recombination l lOsterbrock Il989h . Then the Lyo- luminosity is 
related to the SFR (M SF ) by, 



-'Lya 



0Mhv„(l -f cx )N y M &F . 



(1) 



Here, hv a = 10.2 eV and / csc are, respectively, the energy of 
a Lyman-a photon and the escape fraction of Lyman continuum 
photons. We use f esc ^ 0.1 as meas ured from stacked spectra 
of high-z LBGs ( Sh aplev et all I2006T) . Further, N y is the num- 
ber of ionizing photon s released per baryon of star formation. 
This is mainly a function of th e metallicity and the stellar initial 
mass function (see Table 1 of ISamui et al.l 120071) . As high red- 
shift DL As have generally low m etallicities, [Z/Zq] ~ -1.5 (see 
Fig 3 of iNoterdaeme et alj |2008), we interpolate values of Table 
1 of ISamui et alj l l2007t) to have N y = 9870 corresponding to 
[Z/Z ] = -1.5 and a Salpeter initial mass function with a = 2.35. 
Lya photons are sensitive to different kinds of radiation transport 
effects including resonant scattering (within the ISM as well as the 
intergalactic medium (IGM)) and attenuation via dust. Hence, the 
observed Lyo- luminosity can be related to the emitted luminosity 

by 




i obs 
^Lya 



J esc J -'Lyai 



(2) 
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Figure 2. Evolution of the cosmic star formation ra te density with red- 
shift. T he black points with error bars are taken from Hopk ins & Beacoml 
12006). The red point refers to the 3<x upper limit on the contribution of 
DLA galaxies with log N(H I) > 20.62 as measured in this work. The blue 
point is the same assuming all DLAs with log N(H I) > 20.30 contribute the 
same. The green po ints show the contribution of Mg n systems estimated by 
Mena rd et al] J2009l) . 



where f^." is the escape probability of the Lya photons. It is 
well known that only a small fraction of LBGs are Lyo- emitters 
and even in these ca ses only a small fraction of Lyo- photons es- 
cape the galaxy (see lKornei et al.|[2010h . Taking this fact into ac- 
count it i s appropriate to us e a volumetric / es £ a of 0.05 as esti- 
mated by lHaves et alj j201fj) for high-z galaxies. With these as- 
sumptions, the 3tr upper limit we have obtained above corresponds 
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Figure 3. Wavelength range of the SDSS composite spectrum encompass- 
ing the O vi absorptions. Gaussian fits to some of the absorption lines are 
also overplotted 

to a SFR of < 1.2 M Q yr . For comparison, this limiting value is 
higher than the average SFR of 0.1-0.5 Mm yr -1 measured in low-z 
Mg ii and Ca n sys tems dWild et alJl2007l : lNoterdaeme et al]|201(]| : 
iMenard et all2009h . 

The SFR per unit comoving volume for DLAs is given by 

AJ.Cz) = < >(A/A p )dN/dX (3) 

where < i//»(z) > is the average SFR per unit physical area at z, A 
is the average physical cross-sectional area, A p is the average pro- 
jection of A on the plane of the sky, and dN/dX is the incidence 
of DLAs per unit absorption distance interval. We convert the star 
formation rate (Msf) to < i//,(z) > by assuming the typical size of 
DLAs to be similar to that of faint Lycr emitters, R = 10 kpc. 

We estimate dN/dX to be 9.95 x 10~ s Mpc" 1 and p. < 0.07 
M yr 1 Mpc~ 3 for DLAs with log 2V(Hl) > 20.62 at < z >= 2.86 
assuming (A/A p ) = 2. If we assume the upper limit on the Lya 
flux is also valid for all DLAs with log N(H i)>20.3, then we get 
p, < 0.13 M yr~' Mpc~ 3 . It is to be remembered that increasing 
R or /egj!" will further decrease p,. In Fig.[2]we compare these val- 
ues with the extinction corrected c omoving star formation rat e den- 
sity measurements summarized in Hopkins & Beacom ( 2006). The 
conservative upper limits from our study are comparatively lower 
than the measurements based on LBGs. 



5 HIGH-IONIZATION GAS PROBED BY O vi 
ABSORPTION 

If DLAs are associated with star forming regions, they should be 
associated with hot gas related to galactic outflows and/or infall 
of primordial gas. Indee d, numerical simu lations of cosmological 
structure formation (e.g iDave et al.|[200lh predict that gas falling 
onto diffuse large-scale structures will be accelerated to supersonic 
speeds and become shock-heated to temperatures of 10 5 to 10 7 K, 
creating a phase known as the warm-hot IGM. Similarly a sustained 
star formation activity in galaxies can lead to super novae driven 
outflows of hot gas. Observations of high-z Lyman Break Galaxies 
(LBGs) frequently show galactic scale superwinds and the mass 



outflow rate may b e intimately related to the star formation rate 
jPettini et alj200l|). T his hot gas can be traced by O vi absorption. 

iFox et al. I J2007I) have detected Ovi in at least 34% of the 
DLAs and suggested that Ovi could be present in 100% of them. 
We confirm this finding by detecting the Ovi doublet (with rest 
equivalent width W,(0 iv/11032) = 0.4A) in our composite spec- 
trum (see Fig. [3). The average deconvolved velocity width of the 
Ovi lines (~275 kms~') is systematically higher than that of C iv 
(~250 kms~'), Si iv (~205 kms~') and singly ionization lines 
(~120 kms~'). The velocity range spanned by the hot phase is 
therefore a factor of two larger than that of the neutral phase. This 
suggests a connection between actively star forming regions with 
outflowing and/or infalling hot gas and DLAs. 



6 SUMMARY AND DISCUSSION 

Using a sample of 341 DLAs with log ,/V(H i)>20.62 we place a So- 
upper limit of 3.0xl0~ ls erg s -1 ctrT 2 on the Lya flux emitted at the 
redshift of the DLA within ~1.5 arcsec (or 12 kpc at z~2.86) to the 
QSO line-of-sight. This corresponds to an average hya luminosity 
of < 2x 10 41 erg s or an object with a luminosity of 0.03 L,(Lyci'). 
Thus, the typical DLA population seems to originate from the low 
luminosity end of the high redshift Lyo- emitting galaxies. 

At low redshift the star forming galaxies obey the Kennicutt 
& Schmidt law, 

<iji ir > x = when N ± < Nf' 

= K(NJNl~f„ Nj_ > Nf 1 . (4) 

iKennicuttl Jl998al fbh have found K = (2.5 ± 0.5) x 
10~ 4 MO yr 1 kpc~ 2 ,/? = 1.4 ± 0.15 , Nl = 1.25 x 10 20 cm 2 and 
log N L ™ = 20.62. Using the above equation and the yV(H i) distribu- 
tion of our sample we would expect to measure in our experiment 
a surface SFR of 5.5 x 10~ 3 M Q yr~' kpc -2 . Assuming DLAs have 
a radius of 10 kpc gives an expectation of 1.7 ± 0.4 M Q yr -1 . 
This is higher than our 1.2 M Q yr~' 3cr upper limit derived in the 
previous Section. This suggests that star formation in DLAs is 
less efficient compared to what is seen in the H i disks of nearby 
galaxi es. This is consistent with the conclusions bv lWolfe & Chenl 
(2006). In the local universe such low star formation efficiencies 
are seen either in the outer regions of galaxies or in dwarf galaxies 
dRovchowdhurv et al.ll2009l ; lBigiel et alfeoiOl) 

Using Cn* absorption lines. Iwolfe et aT | d2003al> have es- 
timated the SFR per unit area, = 10~ 219 and 10~ 132 
M Q yr~' kpc~ 2 for CNM (cold neutral medium) and WNM (warm 
neutral medium) models, respectively. The pure WNM solution 
will require an unphysically low escape fraction (< 0.004) of Lycr 
photons in order for the model to comply with our upper limit. As 
DLAs have very little dust content, such low Lya escape fraction is 
unlikely. In the case of the CNM solution, for the upper limit we ob- 
serve to be consistent with the C n* measurements, the Lyo- escape 
fraction / e ^" should be smaller than .03. This value, alt hough low, 
is similar to that estimated in LBGs dHaves et al.ll2010f). Note tha t 
the high-cool population of C n* as defined bv lWolfeetal] ( l2008h 
will also require unphysically low values of /Jc*. In any case, it can 
be concluded that in situ star formation is low in DLAs. 

Using deep long-slit spectroscopy of an empty field, 
iRauch et al.l ( 120081) have detected extended and weak Lyo- emit- 
ters at 2.67< z <3.75 Based on the projected size (a median radius 
of 0.99" corresponding to a physical size of 7.7 kpc) and a vol- 
ume density of 3xlO~ 2 fr 3 Mpc~ 3 , they found that the number per 
unit redshift of these emitters is consistent with that of DLAs. The 
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individual Lyo- flux in these objects range from 1.5 to 15xl0~ 18 
erg s" 1 cm~ 2 with mean and median fluxes of, respectively, 3.7 
and 3xl(T 18 erg s cirT 2 . If typical DLAs were associated with 
these objects (with sizes less than the projected size of the SDSS 
fibre), then we should have detected Lyo- emission in the stacked 
spectrum. Thus the analysis presented here does not support the 
idea that the unresolved or cent rally dominated faint Lyo- emit- 
ters found by iRauch et al.l J2008h are the host-galaxies of DLAs 
with log N(H i) > 20.62. However, our analysis can not rule out 
DLAs being similar to their amorphous extended class of objects 
with sizes > 10 kpc. 

In summary, there are strong evidences for DLAs to be asso- 
ciated with regions of star-formatio n: (i) the median metalli city of 
DLAs is of the order of 10" 15 solar I Noterdaem eet al .120081) and is 
always larger than 10~ 3 solar jPenprase et al.l201Ch : (ii) weconfirm 
the findings by Fox et al. (2007) that strong O vi absorption is as- 
sociated with DLAs with a velocity spread (FWHM ~ 250 km s" 1 ) 
about twice larger than the velocity spread of the low-ionization gas 
suggesting bulk-motions in the hot phase due to kinematical input 
from supernovae and/or gravitational effects; (iii) C n* is detected 
in almost half of the DLAs and can only be explained by energy 
input from star-formation activity (Wolfe et al. 2008). However, we 
show here that the mean Lyo- emission within 12 kpc from DLAs 
is smaller than the emission in most of the weak Lyo- emitters de- 
tected by Rauch et al. (2008) and is consistent with very low in situ 
star formation rates. Therefore, the DLA phase could be spread far 
away from the center of the star forming regions. The hot phase 
associated with DLAs could then be the relics of past wind episods 
in the star forming regions. 
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